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For multilayered materials, reflectivity depends on the complex dielectric function of all the constituent layers, 
and a detailed analysis is required to separate them. Furthermore, for some cases, new quantum states can 
occur at the interface which may change the optical properties of the material. In this paper, we discuss 
various aspects of such analysis, and present a self-consistent iteration procedure, a versatile method to 
extract and separate the complex dielectric function of each individual layer of a multilayered system. As a 
case study, we apply this method to LaAlOs/SrTiOs heterostructure in which we are able to separate the 
effects of the interface from the LaAlC >3 film and the SrTiOs substrate. Our method can be applied to other 
complex multilayered systems with various numbers of layers. 
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I. INTRODUCTION 

Recent technological advances in synthesizing multi¬ 
layered materials with precise atomic control have made 
it possible to study the various exotic quantum phe¬ 
nomena that can occur at the interfaces of dissimilar 
materials 1 . The interplay between charge, spin, and or¬ 
bital at these interfaces can lead to various exciting phe¬ 
nomena such as orbital and spin reconstructions, metal- 
insulator transitions, magneto-electric coupling, super¬ 
conductivity, quantum Hall effect, and topological ef¬ 
fects. What makes it even more interesting is that these 
phenomena can occur even if the parent materials that 
make up the interface are not known to exhibit those 
properties. Examples include the ferromagnetic inter¬ 
face (and associated colossal magnetoresistance) between 
antiferromagnetic insulators LaMnC >3 and SrMnO. 3 2 ’ 3 , 
superconducting interface between insulator La 2 Cu 04 
and metallic overdoped (La,Sr) 2 Cu 04 4 , conducting inter¬ 
face between band insulators LaAlC >3 (LAO) and SrTiOs 
(STO) D , and many others 1 ’ - . Another case is graphene; 
since it usually needs to be suspended on top of a sub¬ 
strate due to its two-dimensionality, it is also effectively 
a multilayered system 8 ’ 19 , and its interaction with the 
underlying substrate can be very interesting to study 19 . 

In order to study the nature and mechanisms behind 
these various interface phenomena, it is very crucial to 
have a thorough understanding of the electronic band 
structure at the interface and how it differs from those 
of the parent materials. One way to directly probe this 
is by measuring the complex dielectric response of the 
material in a broad energy range 18 ’ 20-22 . For example, a 
combination of spectroscopic ellipsometry (0.5 - 5.6 eV) 
and ultraviolet - vacuum ultraviolet (UV-VUV) reflectiv¬ 
ity (3.7 - 35 eV) can be used to obtain the reflectivity of 
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the material in the broad range of 0.5 - 35 eV 18 ’ 20-22 . 
The broad photon energy range is crucial in order to 
yield a stabilized Kramers-Kronig analysis of the reflec¬ 
tivity data (explained in details later), so that the correct 
complex dielectric function e(ui) = £i(w) + *£ 2 ( 0 ;) of the 
material can be extracted reliably from the reflectivity. 
This technique has been proven to be important in the 
study of a wide variety of materials, ranging from man- 
ganites, graphene, and oxides such as SrTiOs 18 ’ 20-23 . 

One important difference between the study of inter¬ 
faces with that of bulk materials is the fact that the in¬ 
terface is buried under one or more layers of parent ma¬ 
terials. Thus, any technique intended to be used in the 
study of interfaces has to be able to probe the buried in¬ 
terface without disturbing the parent materials surround¬ 
ing it. For example, in LaAlOs/SrTiOs (LAO/STO) het¬ 
erostructure, the LaAlOs film thickness is typically in the 
order of 1 - 4 nm, which means that the interface is buried 
also at 1 - 4 nm below the surface. In the optical re¬ 
flectivity setup described above, the photon penetration 
depth is found to be in the order of 10 - 40 nm, which 
is more than sufficient to probe the buried interface of 
LaAlOs/SrTiOs. In general, this is also applicable to 
other multilayered systems as long as the depth at which 
the interface is buried does not exceed the penetration 
depth of the photon. 

In this paper, we discuss various aspects of reflectivity 
and dielectric function analysis of both bulk and multi¬ 
layered materials. Especially for the multilayer analysis, 
we present a self-consistent iteration procedure, a versa¬ 
tile method to extract and separate the complex dielec¬ 
tric function of each individual layer of a multilayered 
system. As the case study, we apply this method to 
LaAlOs/SrTiOs heterostructure ', in which we are able 
to separate the effects of the interface from the LaAlOs 
film and the SrTiOs substrate. Our method can be ap¬ 
plied to other multilayered systems with various numbers 
of layers. 
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II. GENERAL ANALYSIS OF OPTICS DATA OF BULK 
MATERIALS 

We first discuss the spectroscopic ellipsometry (SE), -g 
which can cover a photon energy range of 0.5 - 5.6 eV. SE g 
is a self-normalizing technique to determine the complex aj 
element of dielectric tensor from a single measurement 
without the need to perform a Kramers-Kronig trans¬ 
formation, making it free from any ambiguities that are 
related to the normalization of conventional reflectivity 
results 1 ' ! . The raw data measured SE is expressed in 
terms of T and A, which are defined as 25 

p = tan T exp(i A) = —, (1) 

r s 

where p is the ratio between r p and r s , the reflection 
coefficients of p- (parallel to the plane of incident) and s- 
(perpendicular to the plane of incident) polarized light, 
respectively. From the Fresnel equations, r p and r s are 
defined as 

ij _ rij cos Qj - Tij cosfy 
p rij cos 0i + rii cos 9j 

and 

r ij _ rij cos Qj - rij cos 9j ^ 

s rii cos 0i + rij cos 9j 

Here, n and 9 represent the refraction index and angle 
of incident from the surface normal, respectively. The i 
and j indices represent the two materials involved in the 
photon propagation. 

From n, the complex dielectric function e{uf) can be 
obtained using 

VeR = n(u), (4) 

where u> is the photon frequency. The s{u>) obtained us¬ 
ing Eq. 4 can then be converted back to reflection coef¬ 
ficients using Eqs. 2 and 3. The tanT and A are essen¬ 
tially ratios of the intensities (for tanT) and phases (for 
A) of the reflection coefficient of the p- and s-polarized 
lights (Eq. 1), which makes them (and any quantities de¬ 
rived from them, including reflectivity) self-normalized. 

This is an important attribute of SE, since the converted 
self-normalized reflectivity can be used to normalize the 
reflectivity data obtained using other methods, such as 
the UV-VUV reflectivity technique. 

Next, we discuss the analysis of the UV-VUV reflec¬ 
tivity data (3.7 - 35 eV) of bulk materials. In optics, a 
material can be considered as a bulk material when its 
thickness, d 1 is more than five times the photon penetra¬ 
tion depth, D , (see Section IX for detailed discussion of 
D), i.e. d > 5 D . Due to the self-normalized nature of 
SE, the SE-derived reflectivity can be used to normalize 
the UV-VUV reflectivity at the low energy side within 
the range of 3.7 - 5.6 eV. Furthermore, the high energy 
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FIG. 1. The normalized high-energy reflectivity (0.5 - 35 eV) 
of SrTiOs is compared to the self-normalized reflectivity ob¬ 
tained from spectroscopic ellipsometry (0.5 - 5.6 eV), the un¬ 
normalized UV-VUV reflectivity (3.7 - 35 eV) from the UV- 
VUV reflectivity measurements (scaled down by 20 x to fit 
the graph), and the calculated reflectivity from off-resonance 
considerations (> 30 eV). To obtain the normalized reflectiv¬ 
ity in the full range of 0.5 - 35 eV, the unnormalized UV-VUV 
reflectivity is further scaled down to match the spectroscopic 
ellipsometry and the off-resonance data, and then the three 
data are appended together. The raw data are reproduced 
with permission from T. C. Asmara et al., Nat. Commun. 5, 
3663 (2014) 2| . Copyright 2014 by Nature Publishing Groups 
(NPG). 

part (> 30 cV) is normalized using calculations based on 
off-resonance scattering considerations according to 

r = i^-F(0)P f (29) > (5) 

smti 

where r is the reflection coefficient, rg is the classical elec¬ 
tron radius (e 2 /me 2 ), A is the photon wavelength, Pf(9) 
is the polarization factor (equal to unity for s-polarized 
light and equal to cos 9 for p-polarized light), and F(9) 
is the structure factor per unit area given by 

F(9) = y^n q f q exp(i^flsm9). (6) 

The summation is performed over the different types of 
atoms on a particular atomic plane on which the light is 
incident, with n q denotes the number of atoms of type q 
in that particular plane, f q denotes the tabulated atomic 
form factor corresponding to that atom q, and z q de¬ 
notes the direction vector normal to the plane in ques¬ 
tion. From the above step, normalized reflectivity in the 
range of 0.5 - 35 eV can be obtained. As an illustration 
for the normalization procedure, Figure 1 shows the nor¬ 
malized high-energy reflectivity (0.5 - 35 eV) of SrTiOs 
as compared to the self-normalized reflectivity obtained 
from SE (0.5 - 5.6 eV), the unnormalized UV-VUV reflec¬ 
tivity (3.7 - 35 eV), and the off-resonance considerations 
(> 30 eV). (The details of the experimental procedures 
used to obtain the data are discussed in Section IV.) 

From the above step, normalized reflectivity in the 
range of 0.5 - 35 eV can be obtained. For isotropic 
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bulk materials, the e(u>) can be extracted from the 
normalized UV-VUV reflectivity using Kramers-Kronig 
analysis' 8 ~ . The procedure is as following. The reflec¬ 
tion coefficient and phase difference between the reflected 
and incident light, tp, are related through the Kramers- 
Kronig transformation according to 

r(w) = ^R(uo) exp(ip(uj)) (7) 

and 

tp(u,) = --pf ln R ( x ] ) dx + ip(0), (8) 

7T J 0 X 2 ~ 

where R = |r| 2 is the reflectivity and P is the Cauchy 
principal value. From here, the refractive index n and 
extinction coefficient k can be extracted from reflectivity 
using 


1 + R — 2 y/R cos p 


III. GENERAL OVERVIEW OF ANALYSIS OF 
MULTILAYERED MATERIALS 

If the material is in isotropic bulk form, the Kramers- 
Kronig analysis is straightforward, since the obtained re¬ 
flectivity only depends on e(w) of one material. However 
if the material is composed of several layers (i.e. a mul¬ 
tilayer), the analysis becomes more complex due to the 
interference between the light reflected from the surface 
and those reflected from the interface(s). For example, 
according to analysis of wave propagation in a strati¬ 
fied medium, the reflection coefficient of a thin film on a 
substrate (i.e. a two-layered material) has the following 
form 51 ’, 


^amb, multi 


T‘amb,film “I - T film, subs & 

1 H" T'amb,film^film,subs€}^f tlrn 


(16) 


where 


and 


2 R sin p 

1 + R — 2 sJR cosp 


( 10 ) 


Finally, the real (£i) and imaginary (£ 2 ) parts of complex 
dielectric function can be obtained from n and k via 


and 


2 7 2 

£i = n — k 

(ii) 

£2 = 2 nk. 

(12) 


Like R and p, £\ and £2 are also related through the 
Kramers-Kronig relationship according to 


and 


£i(w) - 1 = -P 

7 r 



X£ 2 (x) 
X 2 — up 


dx 


£2(w) 



£l(x) 47TCT_dc 

—t- ydx H- 

X z — LO z (jJ 


(13) 


(14) 


where ape is the DC conductivity. 

The Kramers-Kronig analysis can be either done di¬ 
rectly through function inversion (i.e. by directly using 
Eqs. 7-12), numerical approximation, or through fitting. 
In this paper, the analysis is done by fitting 35 using 
the Kramers-Kronig-transformable Drude-Lorentz oscil¬ 
lators according to 


£(w) = £00 + 


UJ, 


p,k 


w 0,fc-W 2 -* r few' 


(15) 


The high-frequency dielectric constant is denoted by £ 00 ; 
w p ,fc, u>o,k, and T*, are the plasma frequency, the trans¬ 
verse frequency (eigenfrequency), and the line width 
(scattering rate) of the fc-th oscillator, respectively. Since 
the energy range involved is very broad (covering 0.5 - 
35 eV), the analysis yields a stabilized Kramers-Kronig 
transformation. 


S i = ^n i cos9 i . (17) 

Here, the subscripts amb , multi , /*/m, and subs repre¬ 
sent the ambient, the multilayer, the thin film, and the 
substrate, respectively, which are the various materials 
involved in the propagation of the photon. 

In other words, for a multilayered system the ob¬ 
tained reflection coefficient, reflectivity (R am b,multi = 
\r a mb,muiti\ 2 ), and, via Eq. 1, p, 4>, and A depend on the 
e(ui) of both films and substrate, along with the thickness 
of the film and the angle of incidence 3 ' , i8 , 

Pamb, multi — R(£ film : £ subs , O') , (18) 

Pamb,multi — p(&filmi E-subs: dfUmi ^); 

and a detailed analysis is required to separate them. In 
this paper, we discuss various aspects of such analysis, 
and present the self-consistent iteration procedure, a ver¬ 
satile method to extract and separate the e(w) of each in¬ 
dividual layer of a multilayered system so that they can 
be further analyzed separately. 

If the properties of the substrate are not expected to 
be significantly affected by the presence of the films, the 
reflectivity of the bare substrate can be separately mea¬ 
sured, and from that its e(u) can be separately obtained 
using the general procedure described in Section II. Then, 
if there is only one film layer, the £(w) of the film can 
be straightforwardly obtained from the total reflectivity 
using Eqs. 2, 3, 16, and 17. However, if there are mul¬ 
tiple layers of films composed of different materials, and 
the e(w) of each material is unknown (or different from 
their bulk forms), then the analysis becomes more com¬ 
plicated. This is because there are several unknowns but 
only one equation (Eq. 16) ;8 , which prevents a straight¬ 
forward mathematical solution. The same problem also 
occurs if the properties of the substrate are affected by 
the presence of the films. For example, if parts of the 
substrate near the interface become modified due to the 
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FIG. 2. (a) Crystal structure of LaA 103 /SrTiC> 3 . (b) Trans¬ 
port measurement of LaAlOs/SrTiOa, showing the LaAlC> 3 - 
thickness-dependent metal-insulator transition. The raw data 
are reproduced with permission from T. C. Asmara et al, 
Nat. Commun. 5, 3663 (2014)" . Copyright 2014 by Na¬ 
ture Publishing Groups (NPG). (c) Multilayer consideration 
of LaAlOa/SrTiOs. 

presence of the films, then in the analysis the interface 
needs to be treated as an effective additional layer. So, 
even if the system is composed of only one layer of thin 
film on top of a substrate (an apparent two-layered sys¬ 
tem), it needs to be treated as if it was a three-layered 
one due to the presence of the interface layer. 

IV. CASE STUDY: LAO/STO HETEROSTRUCTURE 

One example is the LaAlOa/SrTiOa heterostructure 
(Figure 2 (a)), which is also the system that will be 
used as the case study throughout this paper. When 
thin film of LaAlOa is deposited on SrTiOa, a conduct¬ 
ing quasi two-dimensional electron gas (2DEG) with high 
mobility and carrier density emerges at the interface . 
Interestingly, this quasi-2DEG only emerges when the 
LaAlOa film reaches a certain critical thickness, usually 
4 unit cells (uc) or more 1 ' 1 . Below the critical thick¬ 
ness, the system remains insulating. Furthermore, the 
interface is also found to exhibit magnetism' 1 ' - ' 1 and 
superconductivity 45- , and two-dimensional coexistence 
of both has even been observed 1 - . These observations 
are very remarkable considering LaAlOa and SrTiOa are 
non-magnetic insulators in their bulk state . Thus, 
the understanding of electronic band structure at the 
LaAlOa/SrTiOa interface is very crucial to reveal the na¬ 
ture and mechanisms of these interesting phenomena. 

For the case study, four samples of LaAlOa/SrTiOa 
with varying thicknesses of LaAlOa (2, 3, 4, and 6 uc, 
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respectively) are prepared using techniques described 
elsewhere" ’ . From transport measurements, it is 
known that the 2 and 3 uc samples are insulating with 
carrier density and conductivity below the measure¬ 
ment limit, while the 4 and 6 uc ones are conducting 
with carrier density of 4 — 6 x 10 13 cm -2 and conduc¬ 
tivity of 4 — 8 x 10 -5 fl -1 , consistent with previous 
results 39,40 ’ 42 ’ 45 (Figure 2 (b)). Bulk LaAlOa and bulk 
SrTiOa samples are also prepared for comparison. 

The reflectivity is obtained using a combination of 
SE (0.5 - 5.6 eV), and UV-VUV reflectivity (3.7 - 35 
eV) ltS ’ 211- . The details of the optical measurements are 
as follow. The SE measurements are performed in the 
spectral range between 0.5 and 5.6 eV by using an SE 850 
ellipsometer at room temperature . Three different inci¬ 
dent angles of 60°, 70°, and 80° from the sample normal 
are needed, and the incident light is 45° linearly polarized 
from the plane of incident. For reflectivity measurements 
in the high-energy range between 3.7 and 35 eV we use 
the superlumi beamline at the DORIS storage ring of 
HASYLAB (DESY) 51 . The incoming photon is incident 
at the angle of 17.5° from the sample normal with linear 
polarization parallel to the sample surface. The sample 
chamber is outfitted with a gold mesh to measure the 
incident photon flux after the slit of the monochromator. 
The measurements are performed in ultrahigh vacuum 
environment (chamber pressure of 10 -9 mbar) at room 
temperature. The obtained UV-VUV reflectivity data is 
calibrated by comparing it with the luminescence yield 
of sodium salicylate (NaCrHsOa) and the gold mesh. 

For the nearly-isotropic bulk LaAlOa and bulk SrTiOa, 
the reflectivity normalization along with the extraction 
of e(cu) from the normalized reflectivity can be performed 
using the general procedure described in Section II. (The 
normalized reflectivity and extracted e(w) of both are 
shown later in Figure 6.) On the other hand, the analy¬ 
sis of LaAlOa/SrTiOa is not as straightforward, due to its 
heterostructure nature as well as the presence of the con¬ 
ducting layer at its interface. For this reason, a multilayer 
consideration based on a boundary analysis of Fresnel 
equation needs to be taken into account in analyzing the 
SE data and the UV-VUV reflectivity of LaAlOa/SrTiOa. 
In this multilayer analysis, LaAlOa/SrTiOa (particularly 
the conducting samples) consists of three layers: the 
LaAlOa film layer at the top, the bulk SrTiOa layer at the 
bottom and an interface layer sandwiched between the 
two at the middle, representing the quasi-2DEG at the 
interface (Figure 2 (c)), consistent with previous obser¬ 
vation using cross-sectional conducting tip atomic force 
microscopy 52 . 

According to analysis of wave propagation in a strat¬ 
ified medium, the reflection coefficient of a three-layer 
system like LaAlOa/SrTiOa can be expressed through 
Fresnel equations as 11 ' 


TambtfLAO T r fLAO,int& f LAO “f* T arnb ^f LAO^fLAO^int^int^STO^ lnt ~\-Tint,STO& ^ (19) 

1 + r amb j L AO r fLAO,inte l2SfLA ° + T f LAO,intfint,STOe- l25ir,t + rambJLAO r int,STO^ l2< ' SfLAO+Sirit ' 1 ’ 
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where the subscripts fLAO , int , and STO represent the 
LaAlOa film, the interface layer, and the SrTiC >3 sub¬ 
strate, respectively. Thus, the reflectivity (and, by exten¬ 
sion via Eq. 1, p, dt, and A) of LaA 103 /SrTiC >3 contains 
mixed information from all three constituent layers 37,38 , 


Ltamb,multi = R{£f LAO, tint, £STO,df LAO i dint ? 
Pamb,multi — pi^ f LA0•> ^inti £ ST0 1 df LAO , dint, 



This makes the extraction of the e(u) of individual 
layer non-trivial, because there are too many unknown 
factors involved. Since e(u) of the bulk SrTiC >3 substrate, 
£sto, can be measured independently and LaA 103 film 
thickness, dfLAO, of each sample are known, Eqs. 1-4, 17, 
and 19 left us with 3 unknown variables: e(ui) of LaAlC >3 
film, £ fLAO, (which might be different from that of bulk 
LaAlOs), e(w) of interface layer, £ int , and the thickness 
of the interface layer, di n t, 

Ltamb, multi — f L AO 5 £ int •> di n t > •> (21) 

Pamb, multi — pij£f LAO, £int, dint, 

(Note that even though e(u>) has real and imaginary 
parts, they are connected through the Kramer-Kronig 
relationship, see Eqs. 13 and 14.) This poses a challenge, 
because (assuming there is no change in £sto across 
the samples) there are 3 unknowns but only 1 equation 
(Eq. 19), which prevents a straightforward mathematical 
solution . To overcome this problem, it can be noted 
that the light phase, S, in Eq. 17 depends mainly on two 
parameters: the incident angle, 0, (angle-dependent)' 53-55 
and the layer thickness, d , (thickness-dependent) 51 . This 
means Eq. 17 can be used to diversify Eq. 19 by varying 
either of these two parameters, so that the number of 
equations can match the number of unknown variables 
(in this case, three). (Note that Eq. 19 can also be di¬ 
versified by varying the ambient within which the mea¬ 
surement is performed 51 ’ 56 , e.g. by purging the measure¬ 
ment chamber with different ambient gas or immersing 
the setup inside different liquids, however concerns about 
surface contamination on surface-sensitive samples may 
make this method less versatile.) This enables us to per¬ 
form a self-consistent iteration procedure on the reflec¬ 
tivity (and thus also p, T, and A via Eq. 1) data, so that 
each unknown variable can be extracted separately. 


V. ANGLE-DEPENDENT ITERATION PROCEDURE 



initial guess: bulk LaAI0 3 
d int is fixed, initial d int = 5 nm 
0 = 60° 


d M is varied 
0 = 70° 


d M is fixed 
0 = 80° 


FIG. 3. Diagram depicting the angle-dependent iteration pro¬ 
cess. 


following (Figure 3). According to previous studies, the 
thickness of the conducting interface might be around 
2-10 nm ,52 ’ 57-61 , so the initial guess for dint, d int \ 0 i, 
can be reasonably set as 5 nm. In general, there are two 
boundary conditions that can be applied when setting the 
initial guess, and also for confirming the physical validity 
of the converged value after iteration, of the thickness of 
films and interfaces: 

1. The thickness should not be lower than the thickness 
of 1 uc of the materials. For films, this is the thinnest 
physical limit for the layer-by-layer deposition, while for 
interfaces this is to take into account any interface rough¬ 
ness effects. 

2. The thickness should not be higher than five times the 
photon penetration depth D (see Section IX for detailed 
discussion of D) 7 since beyond this limit the material is 
optically considered to be bulk- . 

Meanwhile, the initial guess for £jlao can be set as 
the same as e(ui) of bulk LaAlC> 3 , £bLAO , which can be 
obtained independently by measuring it separately. With 
these two variables fixed, Eq. 1 is fitted into experimental 
value of and A measured at 9 = 60° using Eq. 15 
by appropriately adjusting the Drude-Lorentz oscillators 
that make up the e(u) of interface layer ' At this first 
step, by expanding Eq. 21 to the first order, the tentative 
e(ui) of interface after fitting, £i n t[i], can be expressed 


In SE (0.5 - 5.6 eV), the T and A measurements are 
done at three different incident angles: 60°, 70°, and 80° 
from the sample normal, which results in three sets of T 
and A data. Since Eq. 1 can be diversified via Eqs. 17 
and 19 by varying 6 i3-55 , this gives us the three equations 
necessary to perform an angle-dependent iteration pro¬ 
cedure to extract the three unknown variables: £fLAO, 
£int, and d int • 

As a representative, the iteration for the SE data of 
the 4 uc LaAlC> 3 /SrTi 03 sample can be performed as 


£int[ 1] —£int 4“ (dint di n t{ o] ) 

+ (£}LAO — SbLAo) 


& Pamb,multi / 3d 

int 

&Pamb, multi / &£int 
& Pamb,multi/ f LAO 
d pamb,multi/ &£int 


( 22 ) 


It can be seen that at this step, £i n t\i] deviates from 
the actual value of £i n t due to the still-improper values of 
dmt and £fLAO■ To simplify the notation, generalized ad¬ 
dition and subtraction operators, © and 0, respectively, 
can be introduced to represent the correlation effects of 



















6 


dint and e/lao on e int [i} such that, 


£ int[ 1] — £ int © ( £ i n t © £ int{ 0]) © ( £ fLAO © EbLAO )■ (23) 


Here, ef nt and £/ n t[o] are introduced to represent the cor¬ 
relation effects that di n t and d in t [ 0 ] have on e int [ q, re¬ 
spectively. For example, to the first order, £f nt can be 
expressed according to Eq. 22 as, 

d &P a mb : multi/(ddi n t ir,A\ 

£ int ~ di n t • (24j 

VHamb,multi / Uc.%nt 

Another advantage of these generalized operators nota¬ 
tions is that they also allows the higher orders of Eq. 21 
expansion to be implicitly included in Eq. 23, mean¬ 
ing that Eq. 23 can be taken as the generalized form of 
Eq. 22. Thus, due to their convenience, these operators 
shall be used throughout this paper. 

After the first step described above, the newly-fitted 
£ int\ i] is in turn fixed, and di nt is appropriately adjusted 
so that Eq. 1 can now be fitted into experimental value 
of T and A measured at 6 = 70°. Here, the tentative 
value of di n t, dint[i], can be expressed as, 

dint{ i] = di n t © ( d e f LAO 9 d e bLAO ) © ( d\ nt © d e int ^). (25) 

Again, d e fLAO , d \ LAO , d\ nt , and d e int[1] are introduced to 
represent the correlation effects that EfLAO, £ bLAO, £ int , 
and £i n t[ i] have on d int [i ], respectively. As a representa¬ 
tive, d S f LAO can be expressed to the first order as, 


d e fLAO ~ 


dpamb,multi/dEfLAO 
dpamb,multi/&d int 


EfLAO- 


(26) 


Then, the newly-adjusted d,; rlt [l] is also fixed (along 
with the previously-fitted E in t[ i]), and Eq. 1 is fitted into 
experimental value of T and A measured at 9 = 80° us¬ 
ing Eq. 15 by appropriately adjusting the Drude-Lorentz 
oscillators that make up the e(w) of LaAlC >3 film layer. 
Here, the tentative e(oj) of LaAlC >3 film, £/lao[i]> can be 
expressed as, 

e fLAO[l] = EfLAO © (Eint 0 £ int[ 1]) © ( £ int& £ int[l])- ( 2 7 ) 

After that, the process is repeated by going back to and 
A values measured at 9 = 60° and subsequently cycling 
through the incident angles, fitting only one variable at 
each step while keeping the other two fixed. 

Convergence is reached when, at a certain step A, 
£ fLAO[N], £ int[N), and dint[N } can satisfy Eq. 1 for all 
three incident angles, as shown in Figure 4. At this point, 
£ fLAO[N] t EfLAO, £ int[N] t Eint, di n t[N ] ^ dint, and 

the correlation effects between these three parameters are 
minimized (see Section VI for a more rigorous treatment 
of this convergence condition). In other words, the iter¬ 
ation results form a universal fitting that can match the 
data from all incident angles. The iteration thus results 
in the converged values of EfLAO, £ int, and dint, as shown 
in Figure 14 later. Along with the already-known esto 
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FIG. 4. Comparison between the experimentally-measured 
and A of 4 uc LaAlOs/SrTiOs and their fitted values after 
iteration for all three incident angles, (a) For 60° incident 
angle, (b) For 70° incident angle, (c) For 80° incident angle. 
The fitted values match the measured and A very well for 
all three incident angles, confirming the stability of the iter¬ 
ation. The raw data are reproduced with permission from T. 
C. Asmara et al, Nat. Commun. 5, 3663 (2014) . Copyright 

2014 by Nature Publishing Groups (NPG). 


and dfLAO, these quantities can be converted to reflec¬ 
tivity in the 0.5 - 5.6 eV range using Eqs. 2 - 4, 17, and 
19, which then can be used to normalize the UV-VUV 
reflectivity. 

Figure 5 (a) illustrates the iteration process of the 4 uc 
LaAlOs/SrTiOs by showing the evolution of di nt through 
each iteration step. As the iteration progresses, the value 
of dint slowly approaches a distinct asymptotic value, and 
at step 5 it finally converges into ~5.2 nm. 

For 6 uc LaAlOs/SrTiOs, the iteration process can be 
performed similarly, since the only difference is dfLAO, 
which is known and can be appropriately adjusted using 
Eq. 17. Figure 6 shows the fitted 4/ and fitted A after 
convergence that match the measured T and measured 
A of 6 uc LaAlC> 3 /SrTi 03 . The iteration progress of 
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FIG. 5. (a) Iteration progress of the interface layer thick¬ 

ness for 4 uc LaAlOs/SrTiOs. (b) Iteration progress of the 
interface layer thickness for 6 uc LaAlOs/SrTiOs, showing 
the comparison between the arbitrary and predefined starting 
points of 6 nm and 5.2 nm, respectively, (c) Iteration progress 
of the interface layer thickness for 6 uc LaAlOs/SrTiOs, show¬ 
ing the comparison between the arbitrary starting points of 2 
nm and 10 nm, respectively. 


di n t for 6 uc LaAlOa/SrTiOa is shown in Figure 5 (b). 
In this case, the initial guess for dmt is set to be 6 nm, 
and the final converged value is found to be ~5.3 nm, 
very close to the 4 uc value of ~5.2 nm. This indicates 
that the properties of 4 and 6 uc LaAlOs/SrTiOs are 
very similar, and any apparent differences in T, A, and 
reflectivity values between the two samples are mainly 
caused by the difference in dfLAO ■ In fact, because of 
this, since from 4 uc iteration the converged values for 
£fLAO, tint, and di n t of 4 uc LaAlOs/SrTiOs are already 
obtained, those values can also be used as the initial guess 
for the 6 uc LaAlOs/SrTiOs iteration. It can be seen 
from Figure 5 (b) that with those better starting points, 
the iteration process can be simplified and convergence 



FIG. 6. Comparison between the experimentally-measured 'F 
and A of 6 uc LaAlOs/SrTiOs and their fitted values after 
iteration for all three incident angles, (a) For 60° incident 
angle, (b) For 70° incident angle, (c) For 80° incident angle. 
The fitted values match the measured *F and A very well for 
all three incident angles, confirming the stability of the iter¬ 
ation. The raw data are reproduced with permission from T. 
C. Asmara ef al, Nat. Commun. 5, 3663 (2014) . Copyright 
2014 by Nature Publishing Groups (NPG). 


can be achieved with fewer steps, while still reaching the 
same converged value of di nt ~ 5.3 nm. Further tests 
of the stability of the iteration process are also done by 
setting the initial guess for di n t to be 2 nm and 10 nm, 
respectively. From Figure 5 (c), it can be seen that both 
iterations are indeed able to converge to the same di n t 
value of ~5.3 nm, although they need considerably more 
steps to converge because the initial guesses deviate a lot 
more from the converged value. These results confirm 
the self-consistency of the iteration process, showing that 
even if it starts with different initial guesses, the iteration 
does eventually converge into the same final results. 

For the insulating samples (2 and 3 uc 
LaAlOs/SrTiOs), the iteration-based analysis is 
also performed similarly. The fitted Hi and fitted A after 
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FIG. 7. Comparison between the experimentally-measured 'll 
and A of 2 uc LaAlOs/SrTiOs and their fitted values after 
iteration for all three incident angles, (a) For 60° incident 
angle, (b) For 70° incident angle, (c) For 80° incident angle. 
The fitted values match the measured 41 and A very well for 
all three incident angles, confirming the stability of the iter¬ 
ation. The raw data are reproduced with permission from T. 
C. Asmara et al, Nat. Commun. 5, 3663 (2014) . Copyright 
2014 by Nature Publishing Groups (NPG). 



FIG. 8. Comparison between the experimentally-measured 41 
and A of 3 uc LaAlOs/SrTiOs and their fitted values (after 
iteration) for all three incident angles, (a) For 60° incident 
angle, (b) For 70° incident angle, (c) For 80° incident angle. 
The fitted values match the measured 41 and A very well for 
all three incident angles, confirming the stability of the iter¬ 
ation. The raw data are reproduced with permission from T. 
C. Asmara et al, Nat. Commun. 5, 3663 (2014) . Copyright 

2014 by Nature Publishing Groups (NPG). 


convergence that match the measured T and measured 
A of 2 and 3 uc LaAlOs/SrTiOs are shown in Figures 7 
and 8, respectively. For the sake of consistency and to 
make layer-by-layer comparison between insulating and 
conducting LaAlOa/SrTiOa more readily apparent, the 
interface layer is still initially retained in the iteration 
process. However, as shown later in Figure 14, after 
analyzing the normalized reflectivity in the full range of 
0.5 - 35 eV, the e(w) of the (artificial) interface layer is 
found to be very similar to that of bulk SrTiOs, making 
insulating LaAlOs/SrTiOa effectively a two-layer struc¬ 
ture instead. This can be explained by the absence of 
the conducting interface in insulating LaAlOa/SrTiOa. 


VI. THICKNESS-DEPENDENT ITERATION 
PROCEDURE 

From the iteration-based analysis of SE data, the s(oj) 
of each individual constituent layer of LaAlOs/SrTiOa, 
along with their thicknesses, can be extracted. These 
quantities can then be converted into reflectivity within 
the 0.5 - 5.6 eV range using Eqs. 2-4, 17, and 19. From 
here, the normalization procedure of the UV-VUV reflec¬ 
tivity data (3.7 - 35 eV) is similar to that of bulk mate¬ 
rials: using the SE-derived reflectivity to normalize the 
low-energy side (3.7 - 5.6 eV) and the off-resonance scat¬ 
tering considerations according to Eqs. 5 and 6 to normal¬ 
ize the high-energy side (> 30 eV). The normalized re¬ 
flectivity of the LaAlC> 3 /SrTi 03 samples along with that 
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FIG. 9. (a) Reflectivity of LaAlC> 3 /SrTi 03 at different thick¬ 
nesses of LaAlOs film as compared to bulk LaAlOs and bulk 
SrTiOs. (b) Complex dielectric functions, e(to) = £i(w) + 
ie 2 (c>->), of bulk LaAlOs and bulk SrTiOs. The raw data are 
reproduced with permission from T. C. Asmara et al., Nat. 
Commun. 5, 3663 (2014). Copyright 2014 by Nature Pub¬ 
lishing Groups (NPG). 


of bulk LaAlOs and bulk SrTiOs is shown in Figure 9 
(a). 

The same challenge present in the analysis of SE data 
of LaAlOa/SrTiOa due to its multilayered structure is 
also present in analyzing the high-energy reflectivity 
data. Even though dint is already known to be ~5.3 
nm from the SE angle-dependent iteration analysis, it 
still leaves us with two unknowns (high photon energy 
£f lao and high photon energy £*„*), but only one equa¬ 
tion (Eq. 19), which prevents a straight-forward math¬ 
ematical solution. Furthermore, due to a fixed incident 
angle of 17.5° from the sample normal, similar angle- 
dependent iteration as the one done for the SE data can¬ 
not be performed. To address this, we note that Eq. 19 
can also be diversified through Eq. 17 by varying the layer 
thickness 54 , in particular dfLAO • It is for this reason that 
we have intentionally fabricated a pair of insulating sam¬ 
ples (2 and 3 uc of LaAlOs) and a pair of conducting 
samples (4 and 6 uc of LaAlOs). Each pair has similar 
respective physics with only difference in dfLAO, which 
can be rectified by appropriately adjusting Eq. 17. This 
means for each case (insulating and conducting), there 
are two unknowns and two equations for R am b,multi, so 
a self-consistent iteration can be used to extract e(w) of 
each individual layer. 

Let us first discuss the iteration procedure for the in¬ 
sulating samples (2 and 3 uc of LaAlOs/SrTiOs). As the 



FIG. 10. Diagram depicting the iteration process of the 2 and 
3 uc LaAlOs/SrTiOs samples at (a) the initial step (step 0), 
(b) step 1, and (c) step 2. 


initial step (step 0), the 3 uc LaAlOs/SrTiOs is regarded 
as only having two layers: the bulk SrTiOs substrate be¬ 
low and a composite layer (comp) on top (Figure 10 (a)). 
This composite layer represents the mixture between the 
unknown LaAlOs film and the unknown interface layer, 


^comp[0] — EfLAO £int • (28) 

Similar to Eq. 23, Eq. 28 can be explicitly expressed to 
the first order as, 


^comp[ 0] £■ fLAO T Eint 


^ Ramb,multi/ &Eint 
9 Ramb,multi /&£ f LAO 


(29) 


Because the separately-measured esto is known (Fig¬ 
ure 9 (b)), £ comp [o] can be obtained by fitting Eq. 16 into 
reflectivity of 3 uc LaAlOs/SrTiOs (i? 3 UC ) using Eq. 15 
by appropriately adjusting the Drude-Lorentz oscillators 
that make up £ CO mp[ o] (similar to the procedure used in 
the angle-dependent iteration, see Sec. V). After that, 
for the 2 uc LaAlOs/SrTiOs, the interface layer is added 
between the SrTiOs substrate and the composite layer 
so that it now is regarded as having three layers (Fig¬ 
ure 10 (a)). This interface layer is added for the sake of 
consistency and to make layer-by-layer comparison be¬ 
tween insulating and conducting LaAlOs/SrTiOs more 
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readily apparent. The newly-fitted e comp [o] along with 
the already-known £sto is now taken as the input {i.e. 
the composite layer becomes the ” fixed layer”), and then 
Eq. 19 is fitted into reflectivity of 2 uc LaAlOs/SrTiOs 
(R 2 uc) using Eq. 15 to extract the newly-added interface 
layer (the ’’variable layer”) e int r 0 ]. 

For the next step (step 1), the 3 uc sample is also re¬ 
garded as having three layers (Figure 10 (b)), and this 
time the interface layer is designated to be the input fixed 
layer in order to extract e(w) of the now-variable com¬ 
posite layer from 7? 3tlc . After the extraction, the £ CO mp 
for this step 1 becomes 



—» variable layer 
—> fixed layer 

= %ao © o |^ t[0] © 2 & int[j| j 

> 7-1 

£ int[n-l] = £ int[0] © Z 

i =1 

—» fixed layer 
—» variable layer 

£ comp[»] = £ fLAO © £ int ©^ £ int[0] ©Z^ £ int[/] j 




^compfl] £fLAO © £int © £ int[ 0] • (30) 

Then, the focus is again shifted to i? 2 uc- The newly- 
adjusted composite layer (e com p[i]) is designated as the 
input fixed layer to extract e(w) of the variable interface 
layer, 

£ int[ 1] £ int[ 0] © ^mt[ 1] i (HI) 

from i? 2 «c- This £i nt [\] is slightly different than e int [ 0 ], 
by an amount of 5£ int [ 1 ]. In step 2 (Figure 10 (c)), this 
procedure is repeated again, and by the end of that step 
e(ui) of the layers becomes 


£comp[ 2] — £fLAO © £int © i^int[ 0] © ^int[ 1] ) 

(32) 

and 


£int[ 2] £int[ 0] © 1] © ^int[2] • 

(33) 

The iteration procedure is repeated until it eventually 
achieves convergence (see discussion below), and at a cer¬ 
tain general step n (Figure 11), e{ui) of the composite and 
interface layers can be expressed as, 

n —1 

£comp[n\ = £fLAO © £int © i^int[ 0] © ^ ^ ^int[i]) 

i— 1 

(34) 

and 


n 

£int[n\ £int[ 0] © ^ ^ ^^int[i\ m 

i= 1 

(35) 

Eventually, at a certain step N , e(u) of the layers be¬ 
come such that, 

8 £ int{N] 0, 

(36) 

£int[N] ^int[N— 1] ? 

(37) 

and 


^comp[N ] ~ ^comp[N— 1] • 

(38) 

This is our point of convergence. At this point, 
converges to £i n t, 

£int[N] 

N 

£int[N] £int[ 0] © ^ ^ ^ ^int ? 

i=l 

(39) 


FIG. 11. Diagram depicting the iteration process of the 2 and 
3 uc LaAlOa/SrTiOa samples at a certain general step n. 

and e comp [Ar + i] converges to £f L AO , 

N 

£ comp[N- (-1] £fLAO © £ int © (fmt[0] © ^ ^ ^nt[i]) 

i=i (40) 

— t £fLAO © £ int © £ int 
~> £fLAO- 

In other words, when this point is reached, the e(o>) of 
each individual layer is able to be isolated separately, and 
the iteration procedure converges successfully. 

In order to give a better presentation of how this it¬ 
eration procedure is applied to the LaAlOs/SrTiOs case 
study, e(w) of the composite and interface layers of 2 and 
3 uc LaAlC> 3 /SrTi 03 at various iteration steps is shown 
in Figure 12. From there, it can be seen that as the iter¬ 
ation progresses, the difference between e(w) of each con¬ 
secutive step becomes progressively smaller. Eventually, 
the e(u>) at step 6 becomes virtually indistinguishable to 
that of step 5, which means that at step 6 the iteration 
converges. The e(u>) of the interface layer is successfully 
separated from the composite layer, and the e(w) of the 
composite layer becomes equal to the £ (w) of the LaAlC >3 
film layer. 

Moreover, to further ensure the validity of the result¬ 
ing Emt and £fLAO , a consistency check can be per¬ 
formed by inserting e int [jv] and £ corn p[N+i] along with 
£sto into Eq. 19 via Eqs. 2-4, 17, and 19 and confirm¬ 
ing that the resulting R am b,multi can indeed reproduce 
the experimentally-measured reflectivity of both 2 and 3 
uc LaAlOs/SrTiOs simultaneously by appropriately ad¬ 
justing the LaA10 3 thickness factor in Eq. 17 (see Fig¬ 
ure 13). Thus, at the point of convergence the extracted 
e(ui) of each individual layer is able to universally fit the 
reflectivity of both 2 and 3 uc LaAlOs/SrTiOs. 

For the 4 and 6 uc LaA 103 /SrTiC >3 samples (the con¬ 
ducting case), the iteration process should essentially be 
the same as the iteration procedure for the insulating 
case. However, in practice it can be more complex than 
that. This is because in the conducting sample the in¬ 
terface layer is conducting, so its optical properties can 
be very different than that of bulk SrTiC> 3 . This makes 
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FIG. 12. Complex dielectric functions, e(u)) = £\{uj)+i£ 2 {ui), of composite and interface layers of the 2 and 3 uc LaAlOa/SrTiOs 
at various steps as the thickness-dependent iteration progresses, (a) The real part of the dielectric function, ei(w), of the 
composite layer, (b) The imaginary part of the dielectric function, £ 2 ( 0 ;), of the composite layer, (c) The £i(cu) of the interface 
layer, (d) The £ 2 ( 0 ;) of the interface layer. Insets show parts of the plots zoomed at various energy ranges to emphasize the 
evolution of e(w) as the iteration progresses. 


the two-layered structure used in the first part of step 
0 unsuitable to model the strictly three-layered system. 
To circumvent it, this first part of step 0 can be skipped 
entirely. Instead, the 4 uc LaAlOs/SrTiOs is directly re¬ 
garded from the start as having three layers: the compos¬ 
ite layer on top, the interface layer in the middle, and the 
bulk SrTiC >3 substrate at the bottom. The initial guess 
of e(u>) of the composite layer can be set as equal to e(ui) 
of bulk LaAlOa, e C omp[a] = EbLAO- and from here the it¬ 
eration can be continued as normal until convergence is 
achieved and the actual £ 2 ( 0 ;) of both the interface and 
the LaAlC >3 film layers are found. 

Figure 14 shows the e(u) of the composite and interface 
layers of 4 and 6 uc LaAlOs/SrTiOs at various iteration 
steps. Similar with the 2 and 3 uc case, as the iteration 
progresses the difference between e(u>) of each consecu¬ 
tive step becomes progressively smaller, and eventually 
at step 9 the iteration converges and e(cj) of the LaAlOs 
film and the interface layer are successfully separated. 

For consistency check, the fitted reflectivity of the 4 
and 6 uc LaAlOs/SrTiOs after convergence is compared 
to their measured values in Figure 15. It shows that 
the resulting e(oj) of LaAlOs film and interface layer 
are indeed able to closely reproduce the experimentally- 
measured reflectivity of both 4 and 6 uc LaAlOs/SrTiOs 
simultaneously, by appropriately adjusting the LaA 103 
thickness factor in Eq. 17. 


VII. RESULTS: COMPLEX DIELECTRIC FUNCTIONS 
OF LAO/STO 

Figure 16 summarizes the analysis results of LaAlC> 3 - 
thickness-dependent LaAlOs/SrTiOs using the self- 
consistent iteration procedure. It shows the extracted 
e(u>) of LaAlOs film and interface layer for both the 
insulating (2 and 3 uc LaAlOs/SrTiOs) and conduct¬ 
ing (4 and 6 uc LaAlOs/SrTiOs) cases as compared 
to bulk LaAlC >3 and bulk SrTiC> 3 . For LaAlOs film, 
it can be seen (Figures 16 (a) and (b)) that the e(u>) 
of LaA 103 film for both the insulating and conducting 
LaAlOs/SrTiOs are very different than bulk LaAlOs an d 
also as compared to each other. This indicates that there 
are significant differences in band structure and orbital 
occupancy among the different forms (bulk or film) and 
environments (insulating LaAlOs/SrTi 03 or conducting 
LaAlOs/SrTiOs) of LaAlOs, and careful investigation of 
these differences might play a role in revealing the mech¬ 
anisms behind the plethora of interesting phenomena of 

LaA10 3 /SrTi0 3 . 

For the interface layer, Figures 16 (c) and (d) show that 
e(lu) at the interface of insulating samples (2 and 3 uc 
LaAlOs/SrTiOs) is very similar to that of bulk SrTiOs, 
which makes insulating LaAlOs/SrTiOs effectively a two- 
layer system. On the other hand, interestingly for con¬ 
ducting samples (4 and 6 uc LaAlOs/SrTiOs) there is 
a new feature around 8-12 eV for £i(w) and 11 - 16 
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FIG. 13. Fitted reflectivity of insulating LaAlOs/SrTiOs as 
compared to their experimentally-measured values after the 
thickness-dependent iteration convergence, (a) Fitted and 
measured reflectivity of 2 uc LaAlOa/SrTiOa. (b) Fitted and 
measured reflectivity of 3 uc LaAlOa/SrTiOa. The raw data 
are reproduced with permission from T. C. Asmara et al, 
Nat. Commun. 5, 3663 (2014) 2 '. Copyright 2014 by Nature 
Publishing Groups (NPG). 


eV for £ 2 ( 01 ), which can be a characteristic of the 2DEG 
that emerges as the LaAlOs/SrTiOa interface becomes 
conducting. Apart from this new feature, the s(u>) of 
the interface layer of conducting LaA 103 /SrTiC >3 quite 
closely resembles the e(w) of bulk SrTiC> 3 , which indi¬ 
cates that the interface layer is SrTi 03 -like, and that 
the conducting layer mostly resides in SrTiOs side rather 
than LaAlC> 3 . 

To investigate the £(w) of each layer of LaAlOs/SrTiOs 
more thoroughly, knowledge of the band structures of 
LaAlC >3 and SrTiOs is needed to identify the optical tran¬ 
sition of each peak in the e(w) spectra. From this, infor¬ 
mation about the relative changes in orbital occupancy 
and charge transfers among the energy bands that partic¬ 
ipate in the optical transitions can be obtained. However, 
such detailed discussion about the interesting physics of 
LaA 103 /SrTiC >3 is beyond the scope of this paper, and 
thus is covered in our other studies . 


VIII. ERRORS ANALYSIS 


The mean squared errors (MSE) associated with the 
reflectivity fitting process involved in the iteration pro- 


TABLE I. Fitting errors of the thickness-dependent itera¬ 
tion of conducting LaAlOa/SrTiOs. Since R and e are A- 
dependent, \^\fu is a quadratic average over all wavelength 
points, while \5R\fu and \5e\fu are absolute averages over all 
wavelength data points. 


Sample 

1 6-R 1 
\ — \ fit 

\8R\fit 

4 uc LaAlC> 3 /SrTi 03 

4 . 4 % 

0.003 

6 uc LaAlOs/SrTiOs 

9 . 8 % 

0.005 

Average 

7 . 8 % 

0.004 

e 

M fa 


R e[efLAo\ 

0.06 


Im [SfLAo] 

0.09 


Re [Sint] 

0.10 


Im[e int ] 

0.07 



cedure can be estimated using 25 


AR|2 _ 1 \ Rex(Wj) Rfjt(iOj) ,2 

R 'f tt — M — 1 ' Rez(Uj) 

j—l \ jj 


(41) 


where M is the number of measurement points, while 
|^| fa, R ex , and Rfu are the reflectivity fitting errors, 
the experimentally-measured reflectivity, and the fitted 
reflectivity, respectively. To obtain the corresponding 
fitting error for e, \^-\fu, hr the first approximation the 
errors can be propagated using 



it — 


|i 


1 

a/M 


1 ..SR fu 

M 11 R 


(42) 


which is based on Eqs. 2-4. The absolute fitting error 
of e, \Ss\fu is then obtained by multiplying \^-\fa with 
|e| (note that the absolute fitting errors of other quanti¬ 
ties are also estimated this way). For example, Table I 
shows the of the thickness-dependent iteration of 

conducting LaAlOs/SrTiOs and how it propagates into 
\Se\fu. The fitting errors of the thickness-dependent it¬ 
eration of the insulating LaA 103 /SrTiC >3 are similar to 
the values shown in Table I. 

For SE data, since the angle-dependent iteration is per¬ 
formed by fitting the and A data, the corresponding 
|^| f it can be estimated as follows. First, the MSE asoci- 
ated with the fitting of p, |^y|/it, can be estimated based 
on Eq. 1 using 


tip ,2 

7 1 


it 


8 tan Ik 
tan T 



(43) 


where \^^\ 2 fit and |^|^ t are obtained by substitut¬ 
ing the R in Eq. 42 with tan T and A, respectively. Then, 
since p is essentially a ratio of reflection coefficients, in 
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FIG. 14. Complex dielectric functions, e(co) = ei(w)-M£ 2 (w), of composite and interface layers of the 4 and 6 uc LaAlOs/SrTiOs 
at various steps as the thickness-dependent iteration using the modified initial guess progresses, (a) The real part of the dielectric 
function, £i(cu), of the composite layer, (b) The imaginary part of the dielectric function, £ 2 ( 10 ), of the composite layer, (c) 
The £i(cu) of the interface layer, (d) The £ 2 ( 0 ;) of the interface layer. Insets show parts of the plots zoomed at various energy 
ranges to emphasize the evolution of e(oj) as the iteration progresses. 


the first approximation |/it can be approximated to be 
the same as the fitting error for r, \ — From Eq. 43, 

| y \f it can be propagated from |^F| fit using Eq. 42 by 
considering that R = |r | 2 and thus = 2|^| fa- 

For example, Table II shows the of the angle- 

dependent iteration of 6 uc LaAlOa/SrTiOa and how it 
propagates into |fe|/j t . The fitting errors of the angle- 
dependent iteration of other LaAlOa/SrTiOa samples are 
similar to the values shown in Table II. 

Meanwhile, in the first approximation, the fitting error 
of dint , |^ 71 fit, can be propagated from |^| f it using 


, 5 dr 


'I fit — 


\r a mb,LAo\ + \TLAO,int \ + \ fint,STo\ 

21 7"ini,STO |(1 — \ramb,LAo\ 2 ~ \TLAO,int\ 2 ) 

1 ,4b 


14 


'I fit 1 


(44) 


which is based on Eqs. 1 and 19. However, due to the 
thin film limit, \5i n t\ is quite small, which leads to a very 
large \6di n t \fu of several nanometers, comparable to the 
obtained di nt value of ~5.3 nm. On the other hand, Fig¬ 
ures 5 (b) and (c) show that even when several iterations 
are performed with different initial guesses for djnt, they 
are still able to converge accurately to a di nt value of 
~5.3 nm, with a small standard deviation of only ^0.1 
nm. Thus, this means that as long as the initial guess 


TABLE II. Fitting errors of the angle-dependent iteration of 
6 uc LaAlOs/SrTiOs. Since T, A, p, and £ are A-dependent, 

|^|/it, |^ |/it, and |^|fa are quadratic averages over all 
wavelength points, while |54/|/i t , l<5A|/it, and |d£|/it are ab¬ 
solute averages over all wavelength data points. 


Angle 

1 <5^ 1 

1 1 fit 

l**l/« 

1 SA 1 

1 "A" I/it 

l*A|/« 

1^1 fit 

60° 

1.1% 

0.10° 

0.5% 

0.57° 

1.2% 

O 

0 

3.5% 

0.16° 

6.5% 

O 

O 

7.4% 

OO 

0 

0 

1.1% 

0.17° 

8.8% 

0.29° 

8.9% 

Average 

2.2% 

0.14° 

6.3% 

0.65° 

6.7% 

£ 

1^1 m 





Re[£/£Ao] 

0.15 





lm[£ fLAO ] 

0.16 





Re 

0.24 





Im[£i„ t ] 

0.14 






is within a reasonable limit from the actual value of dint 
and the iteration is able to converge successfully, the un¬ 
certainty of di n t obtained from the iteration procedure 
should be much smaller than what \Sdi nt \fit propagated 
from | ^ |/it would otherwise suggest. For this reason, a 
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Photon Energy (eV) 

FIG. 15. Fitted reflectivity of conducting LaAlOs/SrTiOa 
as compared to their experimentally-measured values after 
the thickness-dependent iteration convergence, (a) Fitted and 
measured reflectivity of 4 uc LaAlOa/SrTiOa. (b) Fitted and 
measured reflectivity of 6 uc LaAlOa/SrTiOa. The raw data 
are reproduced with permission from T. C. Asmara et al, 
Nat. Commun. 5, 3663 (2014) 2 '. Copyright 2014 by Nature 
Publishing Groups (NPG). 


good initial guess, for example from the results of other 
techniques, is preferable in order to have a more accurate 
and efficient iteration. This also applies for the uncer¬ 
tainties of other quantities, such as the uncertainties of 
£f lao an d Sint shown in Tables I and II. 


IX. PHOTON PENETRATION DEPTH ANALYSIS 


Photon penetration depth, £), of a material is defined 
as the depth at which the intensity of the incident light, 
/, drops to 1/e of its initial value, Iq , where e is the 
natural constant. It can be obtained from the e(w) of 
the material according to 36 


D= wa. 

2n£2 


(45) 


If the material is multilayered like LaAlOs/SrTiOs, the 
intensity drop depends on the penetration depth of each 
constituent material, in this case the LaA 103 film, the 
interface layer, and the SrTi 03 substrate, 


T / \ T r / dfLAO di n t Z dfLAO di n t 
I{z) = / 0 exp[—(-/-+ —- +--)], 


D f LAO Din 


Dsto 


( 46 ) 


where z > (dfLAO + dmt) is along the direction per¬ 
pendicular to and measured from the surface of the het¬ 
erostructure. From Eq. 46, the effective penetration 
depth, D e ff , of LaAlOs/SrTiOs can thus be expressed 
as, 


Deff =(1 


Dsto 

DfLAO 


)dfLAO + (1 


+ D S to- 


Dsto 

Dint 


) d t 


(47) 


The D e ff of the LaA 103 /SrTiC >3 samples along with 
that of bulk LaAlOs and bulk SrTiC >3 at 17.5° incident 
angle is shown in Figure 17. From the figure, it can be 
seen that the D e ff of all samples is around 10 - 40 nm 
above 5 eV and up to 30 pin below it, which is more than 
sufficient to cover the LaAlC >3 film thickness of 1 - 2 nm 
and the interface thickness of ~ 5 nm. This means that 
the photon can indeed probe the interface thoroughly, 
and even able to penetrate deep into the SrTiC >3 sub¬ 
strate. 


X. EXTENSIONS AND LIMITATIONS OF THE 
ITERATION PROCEDURE 

The iteration procedure is not limited to the analy¬ 
sis of LaA 103 /SrTiC> 3 , and it can be applied to analyze 
various other multilayer systems, even if those systems 
might have different numbers of unknown parameters 
than what is discussed in this paper. For example, if 
three unknown films, each with unknown thickness, are 
deposited on top of a known substrate, there will be six 
unknown parameters in total (three unknown dielectric 
functions and 3 unknown thicknesses). To perform an 
angle-dependent iteration on the system, the SE mea¬ 
surement needs to be done at six different incident an¬ 
gles, so that the number of equations (Eq. 1) can match 
the number of unknowns via Eq. 17. The incident angles 
should be chosen such that the corresponding differences 
in S can give rise to relatively large variations in the mea¬ 
sured ’F and A spectra, so that the iteration can be per¬ 
formed more efficiently. The iteration can then be per¬ 
formed by cycling through the *F and A data measured at 
these 6 incident angles. Of course, with more unknown 
parameters the complexity also increases, which means 
more resources are needed to successfully converge the 
iteration process. 

For the thickness-dependent iteration, assuming that 
the thicknesses of the layers are known from the angle- 
dependent iteration and/or other methods, and can be 
controlled during growth, the above example still leaves 
us with three unknown parameters, which are the s(cu) 
of each of the 3 film layers. To perform the thickness- 
dependent iteration, we need to also prepare three sam¬ 
ples with slightly different thicknesses of those 3 films. 
Again, to make the iteration procedure more efficient, 
the thicknesses should be chosen such that the corre¬ 
sponding differences in 6 can give rise to relatively large 
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FIG. 16. Complex dielectric functions, b(lj) = £i(u) + i£ 2 (ca), of each layer of 2, 3, 4, and 6 uc LaA 103 /SrTiC> 3 , extracted 
from reflectivity using the self-consistent iteration procedure, (a) The real part, £i(<u), of dielectric function of LaAlOs Him, 
as compared to bulk LaA10s3. (b) The imaginary part, £ 2 ( 0 ;), of dielectric function of LaAlOs film, as compared to bulk 
LaAlOa. (c) The £i(a;) of the interface layer, as compared to bulk SrTiOs. (d) The £ 2 ( 0 ;) of the interface layer, as compared 
to bulk SrTiC> 3 . Note that the e(u>) of the LaAlOs hhn of the 2 and 3 uc LaAlOs/SrTiOs are identical due to the nature of the 
thickness-dependent iteration process. The same is true for the £(tu) of the interface layer and for the 4 and 6 uc LaAlOs/SrTiOs 
case. The raw data are reproduced from Asmara et al: . 



FIG. 17. Effective photon penetration depth of 2, 3, 4, and 6 
uc LaAlOs/SrTiOs at 17.5° incident angle along with that of 
bulk LaAlOs and bulk SrTiOs. 


variations in the measured reflectivity spectra. The it¬ 
eration can then be performed by cycling through the 
reflectivity of those three samples, because the number 
of equations (equivalent of Eqs. 16 and 19, extended to 
four layers) now matches the number of unknown param¬ 
eters via Eq. 17. 

However, despite its potentials, there are still some lim¬ 


itations inherent especially to the thickness-dependent 
iteration method. The first is sample variance. Since 
multiple samples are needed to perform the thickness- 
dependent iteration, slight differences in properties 
among the samples (for example due to slightly different 
growth conditions each time the samples are prepared) 
can accumulate to rather large deviations. Because of 
this, extra care is needed to ensure that each sample 
is prepared within almost identical environment. Sec¬ 
ondly, it also needs to be ensured that any variation in 
reflectivity among the samples is only due to the differ¬ 
ent film thicknesses involved ( i.e . only due to Eq. 17). 
Film thickness differences should not significantly modi¬ 
fied the internal properties of the samples, because other¬ 
wise it will render the thickness-dependent iteration pro¬ 
cedure invalid. This is why the iteration cannot be per¬ 
formed between the insulating 3 uc LaAlOs/SrTiOs and 
the conducting 4 uc LaAlOs/SrTiOs, since their inher¬ 
ent properties are modified by the increase of the LaAlOs 
thickness. Furthermore, due to this requirement, the re¬ 
sulting e(w) of each layer is identical for all the sam¬ 
ples involved in a particular thickness-dependent itera¬ 
tion. For instance, e(w) of the LaAlOs film of the 2 and 
3 uc LaAlOs/SrTiOs in Figure 16 are identical to each 
other. The same is true for the e(u>) of the interface layer 
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and for the 4 and 6 uc LaAlOa/SrTiOa case. 

These two limitations are not applicable to the angle- 
dependent iteration, because the measurements at the 
different incident angles are still performed on the same 
sample each time. This eliminates the concern about 
sample variance, and because the iteration is performed 
on data measured from only one sample, the converged 
results are also unique to that particular sample. 


XI. CONCLUSIONS 

In conclusion, we have discussed various aspects of op¬ 
tical analysis of bulk and multilayered materials, as sum¬ 
marized by the flowchart in Figure 18. For analysis of 
multilayered materials, we present the self-consistent it¬ 
eration procedure as a useful tool to separate and extract 
the dielectric functions of each individual layer. The 
method is based on the two main variables that affect 
the reflectivity of a multilayered system: photon incident 
angle (angle-dependent) and layer thickness (thickness- 
dependent). By measuring the samples at different inci¬ 
dent angles or on samples with slightly different layer 
thicknesses, self-converged iteration can be performed 
by cycling through these differently-measured spectra. 
With enough iteration steps, stabilized separation and 
extraction of dielectric function of each individual layer 
can be achieved. By applying the procedure into spec¬ 
troscopic ellipsometry and UV-VUV reflectivity data of 
LaA 103 /SrTiC> 3 , we are able to separate the effects of 
the interface layer from the LaA 103 film and the SrTiC >3 
substrate. With proper adjustments, this method can 
be extended to other multilayered material systems with 
various numbers of layers, making it a very versatile tool 
in analyzing the optical properties of various multilayered 
systems. 
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